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Summary. I discuss the dynamical interaction of galactic disks with the surround- 
ing dark matter halos. In particular it is demonstrated that if the self-gravitating 
shearing sheet, a model of a patch of a galactic disk, is embedded in a live dark halo, 
this has a strong effect on the dynamics of density waves in the sheet. I describe 
how the density waves and the halo interact via halo particles either on orbits in 
resonance with the wave or on non-resonant orbits. Contrary to expectation the 
presence of the halo leads to a very considerable enhancement of the amplitudes of 
the density waves in the shearing sheet. This effect appears to be the equivalent of 
the recently reported enhanced growth of bars in numerically simulated stellar disks 
embedded in live dark halos. Finally I discuss the counterparts of the perturbations 
of the disk in the dark halo. 



1 Introduction 

Dark halos are usually thought to stabilize galactic disks against non-axisymmetric 
instabilities. This was first proposed by Ostriker & Peebles (1973) on the basis 
of - low-resolution - numerical simulations. Their physical argument was that 
the presence of a dark halo reduces the destabilizing self-gravity of the disks. 
Doubts about an entirely passive role of dark halos were raised by Toomre 
(1977), but he (Toomre 1981) also pointed out that a dense core of a dark halo 
may cut the feed-back loop of the corotation amplifier of bars or spiral den- 
sity waves and suppress thus their growth. Recent high-resolution numerical 
simulations by Athanassoula (2002, 2003), also inherent in the work of Debat- 
tista & Sellwood (2000), have shown that quite the reverse, a destabilization 
of disks immersed in dark halos, might be actually true. Athanassoula (2002) 
demonstrated clearly that much stronger bars grow in the simulations if the 
disk is embedded in a live dark halo instead of a static halo potential. This is 
attributed to angular momentum transfer from the bar to the halo via halo 
particles on resonant orbits. Angular momentum exchange between disk and 
halo has been addressed since the pioneering work of Weinberg (1985) in many 
studies theoretically or by numerical simulations and I refer to Athanassoula 
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(2003) for an overview of the literature. Toomre (1981) has shown how the 
bar instability can be understood as an interference of spiral density waves 
in a resonance cavity between the eorotation amplifier and an inner reflector 
(cf. also Fuchs 2004b). Thus it is to be expected that a live dark halo will be 
also responsive to spiral density waves and develop wakes. That this is indeed 
the case has been demonstrated by Fuehs (2004a) employing the shearing 
sheet model. This adds to the confidence in the results of the numerical work 
on bar growth in galactic disks. 

The shearing sheet (Goldreich & Lynden-Bell 1965, Julian & Toomre 1966) 
model has been developed as a tool to study the dynamics of galactic disks and 
is particularly well suited to describe theoretically the dynamical mechanisms 
responsible for the formation of spiral arms. For the sake of simplicity the 
model describes only the dynamics of a patch of a galactic disk. It is assumed 
to be infinitesimally thin and its radial size is assumed to be much smaller 
than the disk. Polar coordinates can be therefore rectified to pseudo-Cartesian 
coordinates and the velocity field of the differential rotation of the disk can be 
approximated by a linear shear flow. These simplifications allow an analytical 
treatment of the problem, which helps also in the present case to clarify the 
underlying physical processes operating in the disk. 

2 Sheeiring sheet model 

The basic disk model is the stellardynamical shearing sheet, which describes 
the local dynamics of a patch of a thin, differentially rotating stellar disk. 
Stellar orbits are calculated in a frame at a distance ro from the galactic 
center, rotating with an angular velocity l2o- Pseudo-Cartesian coordinates x 
and y point in the radial direction and tangential to the direction of galactic 
rotation, respectively. The differential rotation of the disk is approximated as 
a parallel shear flow, v = —2Ax, with A denoting Oort's constant. The surface 
density Eq is assumed to be constant over the entire region. As is well known 
(cf. Julian & Toomre 1966) the stellar orbits in this model are simply epicyclic 
orbits and the phase space distribution function of the stars /o, as derived from 
the time-independent Boltzmaini-equation, is a Schwarzschild-distribution. A 
cartoon of the shearing sheet model is shown in Fig. 1. 



and the response of the disk, /i, is calculated from the linearized Boltzmann- 
equation 



Fig. 1. The shearing sheet model 



The disk is subjected to potential perturbations 




(1) 
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^ + lfo,S<P] + lh^Ho]=0, (2) 

written in general form with Poisson brackets. Hq denotes the Hamilto- 
nian of the stellar orbits in the unperturbed disk. In order to obtain self- 
consistent perturbations the response density has to be inserted into the 
Poisson-equation 

A5^ = 4nG j (fvh, (3) 

where G is the constant of gravitation. Unfortunately, the disk response to 
a single Fourier component of the potential perturbation (1) is not a Fourier 
component of the general disk response. I follow therefore Kalnajs (1971) and 
take scalar products, 

^ j dx j (4) 

of both sides of the Poisson-equation with conjugate basis functions of the 
Fourier transform. Details of the evaluation of the quadratures, which are 
carried out using action and angle variables are given in (Fuchs 2001). The 
result is that the Poisson-equation is converted to an integral equation of 
Volterra-type, which is equivalent to the integral equation given by Julian & 
Toomre (1966), although it is not formulated in shearing coordinates, 

«^k',a; = / dk^K{k^,k'^,k'^,uj)<^\^^y^^uj, (5) 

— oo 

with a kernel /C that can be expressed analytically (Fuchs 2001)^. By Pred- 
holm discretization equation (5) can be transformed into a set of algebraic 
equations. The kernel K vanishes on the diagonal, so that the triangular co- 
efficient matrix of this set of equations has an unity diagonal, implying a 
non- vanishing determinant. Thus the homogenous integral equation (5) has 
no eigensolutions, indicating that there exist in the shearing sheet no - except 
ringlike (fcy = 0) - proper spiral modes in the sense of rigidly rotating spatial 
patterns with well defined growth rates. External potential perturbations or 
initial {t = 0) density or velocity perturbations of the basic state of the disk 
are represented by an inhomogeneous term in the integral equation (5), or the 
corresponding integral equation for the surface density 

^k'.a; = / *x/C(fcx,A;x,fcy,w)Z'k^,k;,aj + rk',a;- (6) 
J — oo 

The resolvent kernel T^ik^, k'^, k'y,u) of the inhomogeneous integral equation 

(6) can be obtained as a Neumann series. Solutions of equation (6) are then 
found in a unique way by a convolution of TZ with the inhomogeneous terms 



^ Positive wave numbers k' will be assumed in the following. 
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rk' • Transforming these solutfons back from io- to time-domain one can show 

that the resulting spatial pattern is a superposition of 'swinging' density 
waves, shearing with the general flow, but exhibiting transient growth as they 
swing by until they finally decay. Fig. 2, taken from Fuchs (1991), illustrates 
this for single plane sinusoidal waves, rt' oc 5{k'y^ — fc™), all of the same initial 
radial wave number A;™. The resulting spatial pattern evolves then as 



Fig. 2. Amplitudes of swing amplified density waves with initial radial wave numbers 
/ci" = — 1 (leading waves). Wave numbers are given in units of fecrit = K^/27rG17o, 

where K denotes the epicyclic frequency. Time is in units of {2Akyy^. The Toomre 
stability parameter is Q = 1.4, and an Oort constant ol A = Qq/2 \a assumed. 



S{x, y, t) = S{t)e'^''^"''+Ky) + n{k^, + 2Afc;t)e'[('=^"+2^'=^*)^+'=>l , (7) 

where the swinging around of the wave crests of the density waves is described 
by the growth of the effective radial wave number fc™ + 2j4fcyt with time. As is 
well known, amplification is high for density waves, which are initially leading, 
but low for initially trailing waves. 

A visual impression of density waves growing from initial random fluctua- 
tions of the surface density of the sheet is given in Fig 3., where snapshots of 
a numerical simulation of the dynamical evolution of the shearing sheet are 
shown (Fuchs, Dettbarn & Tsuchiya, in preparation). 



Fig. 3. Snapshots of the growth of density waves in a numerical simulation of 
the dynamical evolution of the shearing sheet. The sheet was seeded initially with 
random noise as shown in the left panel. The right panel shows the sheet after 
an elapsed time of one epicyclic period. Length unit is the critical wave length 

Acrit = 2-n/kcTii,- 



3 The shearing sheet immersed in a hve dark halo 

The evolution of the distribution function of the disk stars in phase space is 
described by the linearized Boltzmann equation 

dt dx dy dx du dy dv 

d^di + ^hi dUo d^di + ^hi dfdo 



dx du dy dv 



= 0, (8) 
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where {u, v) are the velocity components corresponding to the x and y co- 
ordinates, respectively. Equation (8) has been derived from the general 6- 
dimensional Boltzmann equation assuming delta-function like dependencies 
of the distribution function on the vertical z coordinate and the vertical w ve- 
locity component, respectively, and integrating the Boltzmann equation with 
respect to them. A perturbation Ansatz of the form 

fd = /do + /di , ^d = ^do + ^di , = ^ho + ^hi (9) 

is chosen for the distribution function of the disk stars and the gravitational 
potentials of the disk and the halo, respectively, and the Boltzmann equation 
(8) has been linearized accordingly. 

Similarly the linearized Boltzmann equation for the halo particles can be 
written as 

-|- U— h V— h W- 



dt dx dy dz dx du dy dv 

d^do + '^ho dfhi _ d(I\i + <I\n dfho _ d<Phi + '^di ^/ho 
dz dw dx du dy dv 

a^hi + ^di 5/ho 



dz dw 



0. (10) 



The choice of the dark halo model was lead by the following considerations. 
One of the dc;epcr reasons for the success of the infinite shearing sheet model 
to describe spiral density waves realistically is the rapid convergence of the 
Poisson integral in self-gravitating disks (Julian & Toomre 1966). Consider, 
for example, the potential of a sinusoidal density perturbation 

^ y_oo i-oo ^{x-x'Y + {y-y'f' ^ ' 

^„ /, N , Si(A;a;L) 27rGZ'io sin (fca;) , , 
mx, y = 0)= -4GEio sin (kx) lim \ ^ = ^ — ^ . (12) 

XL^OO k k 

The sine integral in equation (12) converges so rapidly that it reaches at 
kxh = f already 87% of its asymptotic value. Thus the "effective range" of 
gravity is about only a quarter of a wave length. The shearing sheet models 
effectively patches of galactic disks of such size. The wave lengths of density 
waves are of the order of the critical wave length 

- IT" - — Z2 — ' (1-^) 

"-crit rv 

where k denotes the epicyclic frequency of the stellar orbits and Sq the surface 
density of the disk. In the solar neighbourhood in the Milky Way, for instance, 
the critical wave length is Acrit = 5 kpc. Thus it is reasonable to neglect over 
such length scales, like in the shearing sheet model, the curvature of the mean 
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circular orbits of the stars around the galactic center or the gradient of the 
surface density. The curvature of the stellar orbits due to the epicyclic motions 
of the stars, on the other hand, cannot be neglected and is indeed not neglected 
in the shearing sheet model. The radial size of an epicycle is approximately 
given by (Tu/k, where denotes the radial velocity dispersion of the stars, 
and the ratio of epicycle size and critical wave length is given by 

= 0.085Q (14) 

in terms of the Toomre (1964) stability parameter Q which is typically of the 
order of 1 to 2. Concurrent to these approximations I have assumed a dark 
halo which is homogeneous in its unperturbed state. Accordingly the curva- 
ture of the unperturbed orbits of the halo particles is neglected on the scales 
considered here and the particles are assmued to be on straight-line orbits. 
The equations of motion of the halo particles are the characteristics of the 
Boltzmann equation. In a homogeneous halo x = V(^do + ^^ho) = 0, and in 
accordance with this assumption I neglect the force terms V^do and V^ho 
in the Boltzmann equation (10). This simplifies its solution considerably. The 
disadvantage of such a model is that there are no higher-order resonances 
of the orbits of the halo particles with the density waves as described by 
Weinberg (1985) or observed in the high-resolution simulations by Athanas- 
soula (2002, 2003). However their effect was shown to be much less important 
than the main resonances of the particles with the density waves, which are 
properly described in the present model. 



Fig. 4. Sketch of the disk and halo model. 



3.1 Halo dynamics 

The Boltzmann equation (10) can be viewed as a linear partial differential 
equation for the perturbation of the distribution function of the halo particles, 
/hi. with inhomogcneitics depending on the perturbations of the gravitational 
potentials of the disk and the halo, V<?di and V<?hi, respectively. Thus the 
equation can be solved for the disk and halo inhomogeneities separately and 
the solutions combined afterwards by superposition. It was shown by Fuchs 
(2004a) that the Boltzmann equation with the inhomogeneity V??hi describes 
just the Jeans instability of the dark halo. However, dark halos are thought to 
be dynamically hot systems and their Jeans lengths will be of the order of the 
size of the halos themselves. Thus this part of the solution of the Boltzmann 
equation (10) is uninteresting in the present context and will be not considered 
in the following. 

More interesting is the remaining part of the Boltzmann equation (10), 
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+ + + (15) 

dt dx dy dz 



dx du dy dv dz dw 



= 0, 



which describes the halo response to a perturbation in the disk. If the grav- 
itational potential perturbation of the disk is Fourier expanded the Fourier 
terms have the form (cf. equation 33 of Fuchs 2001) 

^dk|,e*('^*+'''=^+''^^)"''iil''l (16) 



with fc|| = |k||| = Y'fcx + fcy. This can be converted to Fourier coefficients of 
the halo potential in 3-dimensional k-space as 

1 r°° 1 fell 

Notice that the coordinate y, which is defined in the reference system of the 
disk, is related to the y coordinate in the reference system of the halo due to 
the motion of the center of the shearing sheet as 

y^y-ro^ot- (18) 

The further solution of the Boltzmann equation is straightforward and is de- 
scribed in full detail in Fuchs (2004a). The distribution function is then inte- 
grated over velocity space to obtain the Fourier coefficients phk of the density 
perturbation of the dark halo. Next the gravitational potential associated with 
this density distribution is calculated from the Poisson equation, 

- k'^'P^k = 47rGphk . (19) 

Since the gravitational forces in equation (8) have to be taken at the midplane 
= 0, it is necessary to convert the solution of the Boltzmann equation ^hk 
from the representation in k space to a mixed representation in (k| | , z) space 
leading to two contributions 

^HL„(-0) = £.fc.^^^<P..l^ (20) 
x{l + tVn ^^^^ erf^* ^^^^ J exp ^^2^, | 



due to halo particles not in resonance with the potential perturbation and 



^Lk,.(-o) = -£*^(^^^^-^v^- 

^ w-kyroOo {kyro[2o - uj)^ 
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due to non-resonant halo particles, pb and cth denote the mean spatial density 
and the velocity dispersion of the halo particles, respectively. The final result 
can be formally written as 

^hk, I = 0) = r (w - fcy ro r^o , fc| I )#dk| I , (22) 

where the real and imaginary parts of T are defined by equations (20) and 
(21), respectively. Thus for any given frequency there is a contribution both 
from the non-resonant and the resonant halo particles. 



3.2 Disk dynamics 

The halo response (22) to the perturbation in the disk has to be inserted into 

equation (8). Solving the Boltzmann equation (8) is greatly facilitated by the 
fact that its form is identical to the case of an isolated shearing sheet with 
the replacement 

^dk ^ (1 + r)^dk (23) 

and one can use directly the results of Fuchs (2001) even though the Boltz- 
mann equation is treated there using action and angle variables instead of 
the Cartesian coordinates as in equation (8). In particular the factor 1 + T is 
carried straightforward through to the fundamental Volterra integral equation 
(equation 68 of Fuchs 2001) 

fK 

^k',c.= / dA;x/C(fcx,A;x)(l+5"(fcx,fc;,a;))^k,,k;,c.+rk',c., (24) 

where the kernel K, is given by equation (67) of Fuchs (2001). rk' describes 
an inhomogeneity of equation (24) related to an initial non-equilibrium state 
of the shearing sheet. Equation (24) is separating in the circumferential wave 
number fcy. In equation (24) the wave numbers are expressed in units of the 
critical wave number fccrit • This implies that the Volterra equation describing 
a shearing sheet embedded in a rigid halo potential is formally the same as 
that of an isolated shearing sheet, because in this case T = and the halo 
mass affects only the numerical values of the critical wave number fccrit and 
the stability parameter Q. It is advantageous to consider equation (24) trans- 
formed back from frequency to time domain. Splitting off the lo dependent 
term expio; 2Ak'' from the kernel and making use of the convolution theorem 
of the Fourier transform of products of two functions leads to 



/-oo / k — k' 



2Ak'y 




(25) 
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where the operator T denotes the Fourier transform from a; to time domain. 
In equation (33) I have assumed an initial perturbation of the disk at time 
t = so that ?^kx,k^,t'<o = 0. The Fourier transform T is given by equation 
(34) of Fuchs (2004a). If this is inserted into equation (25) it takes the form 



^k',t = / " dkji (fcx, fc^) , , ^^..-.i (26) 

+ '^'^ dt'#kx,k^,t'^ |^r(fcx, fc;,a;)e"^^ j } + rk',t • 

Equation (26) can be integrated numerically with very modest numerical ef- 
fort. In Fig. 5 I illustrate the response of the shearing sheet now embedded 
in a live halo to an initial sinusoidal perturbation of unit amplitude. For this 
purpose I use the inhomogeneity term of the Volterra equation 



r-k'. 



/ dA;x£(A;x,A;4)/kx,kv(0) (27) 

J — oo 



derived in Fuchs (2001) with /kx,k'y(0) ^ ^{^^ — k™). The response of the 
shearing sheet to this initial impulse is a swing amplification event as described 
in section (2). The radial wave number fcx evolves as 

fcx = k^' + 2Ak'yt , (28) 

while the circumferential wave number fcy is constant, which means that the 
wave crests swing around from leading to trailing orientation during the ampli- 
fication phase. Around t = 6 the amplitudes become negative which indicates 
that the swing amplified density wave is also oscillating. As can be seen from 
Fig. 5 comparing the evolution of shearing sheets embedded either in a rigid 
halo potential or a live halo this characteristic behaviour of the density wave 
is not changed by the responsive halo, but the maximiuxi growth factor of the 
amplitude of the wave is enhanced by a surprisingly large amount. 



Fig. 5. Swing amplified density wave in the shearing sheet. The upper diagram 
shows the evolution in a shearing sheet embedded in a static halo potential triggered 
by an impulse with unit amplitude and wave vector k'" = (— 2, 0.5)fccrit. Time is 
given in units of (2Afcy°/fccrit)~^. The middle diagram shows the evolution of a 
shearing sheet embedded in a live dark halo triggered by the same impulse. The 
lower diagram shows the difference. The model parameters are A/f^o = 0.5, Q = 
1.4, cTd : iTh = 1 : 5, Gp^/K^ = 0.01, and rofio : CTd = 220 : 44. 



The enhanced maximum growth factor of swing amplified density waves 
due to a responsive halo seems to be the equivalent of the enhanced growth of 
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bars of stellar disks embedded in live dark halos seen in the numerical simula- 
tions. However, the interaction of the shearing sheet and the surrounding halo 
is not only mediated by the resonant halo particles, but the non-resonant halo 
particles play an important role as well. The amplification of density waves de- 
pends critically on the Toomrc stability parameter Q. This is illustrated in 
Fig. 6 where the response of the shearing sheet to the same initial impulse 
as in the previous example is shown, but assuming a stability parameter of 
Q = 2. As can be seen in Fig. 6 there is neither effective amplification of 
density waves in a shearing sheet in a rigid halo potential or in a shearing 
sheet embedded in a live dark halo. 



Fig. 6. Same as in Fig. 5, but adopting Q = 2. 



4 Wakes in dark matter halos 

The perturbations of the gravitational potential and the stirface density of the 
shearing sheet have their counterparts in the dark matter halo. Prom equation 
(19) one obtains 

dK^^<^\.]>.e'''-\ (29) 
where the Fourier coefficients #hk derived in section (3.1) have to be inserted, 
Phk„(.,c.) = -|_^dfc.^i^^>dk||^e^'=- (30) 

Equation (30) can be Fourier transformed back from frequency to time domain 
in the same way as equation (32) of Fuchs (2004a) was Fourier transformed 
back to time domain by a convolution of ^dk||(t') with the Fourier transform 
of the remaining terms under the integral with respect to k^. To these terms 
the operator 

/■oo 

^gi(a,-fc^rof2o)* (31) 



f 

J — ( 



is applied and the integral over w is evaluated with the help of formulae (6.317) 
and (3.952) of Gradshteyn & Ryzhik (2000) leading effectively to expression 

(34) of Fuchs (2004a) with t-t'+^^ replaced by t-t' only. The integral over 

kz can be then calculated using formulae (3.723) and (3.896) of Gradshteyn 
& Ryzhik (2000) with the result 
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Phkll(2,i) = -^{27r^>dk||(<)e"'"''' (32) 
rt _ 

J —oo 

and finally I convert the potential perturbation of the disk to the perturbation 
of its surface density with the relation A;||^dk|| = ~27rGX'dk|| (cf- Fuchs 2001) 
and obtain 




The time integral in equation (33) has been calculated numerically. At the 
time of maximal amplification of the density wave I find, adopting the same 
parameters used to calculate Fig. 5, Phk/Pblpcak = 0.6 Tdk/^olpeafe- Fig. 7 
shows the vertical profile of the density perturbation of the dark halo according 
to equation (33). As can be seen from Fig. 7 the profile indicates a density 
enhancement of dark halo matter close above the density wave crest. At larger 
distances above the midplane there is a density deficit above the density wave 
crest, which means that the cloud of dark matter particles is lagging behind 
like a trail of smoke. 




0.5 1 1.5 



Fig. 7. Vertical profile of the density perturbation in the dark halo induced by 

a density wave in the sheet, which is at maximal amplification, exactly above the 
density wave crest. The relative density contrast Phk/Pb|peak is given in terms of the 
relative density contrast of the surface density of the sheet £'dk/.S'o|peak. 
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The existence of wakes in dark matter halos might have quite practical 
implications. For instance the bar in the Milky Way focusscs dark matter 
particles dynamically into some regions in phase space and depopulates others, 
which may very well affect the flux of dark matter particles through detectors 
on Earth. Such features have been observed as 'star streams' in velocity space 
among the stars in the solar neighbourhood (Dehnen 2000, Miihlbauer & 
Dehnen 2003). These stars are, of course, stars members of the Milky Way 
disk, but similar effects are to be expected among halo objects. 
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